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Abstract 

^^ \ The elements of the quantum mechanical diffusion matrix, leading to a Gibbs equilibrium state 

J>-,' for a set of A^ coupled quantum harmonic oscillators are derived within Lindblad's axiomatic 

approach. Consequences of the fundamental constraints on the quantum friction coefficients are 

cn ' discussed. We derive the equations of motion for the expectation values and variances, and we 

solve them analytically. We apply our results to the description of the charge and mass asymmetry 

p ■ coordinates in heavy- ion collisions, and we investigate the effect of dissipation on tunneling in 

rj ! sub-barrier processes. 

O^. PACS numbers: 03.65.-w; 05.30.-d; 24.60.-k; 25.70.-z 
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I. INTRODUCTION 

Over the last decades, there has been an increasing interest in dissipative phenomena 
in heavy-ion colhsions and nuclear fission [l]. The description of these processes may be 
achieved by the characterization of the dynamics of certain collective degrees of freedom 
such as the relative motion, shape deformation, and mass and charge asymmetries 2|. The 
energy dissipation occurs mainly as a result of the coupling of the slow collective degrees of 
freedom to the fast intrinsic degrees of freedom (e.g. nucleons motion). The latter may be 
regarded as a heat bath in thermal equilibrium at each stage of the reaction . 

Classically, the process of dissipation is quite well understood; a typical example is the 
Brownian motion which has attracted much attention since Eisentien's seminal work. The 
conventional way of treating dissipation and damping of classical systems may be fulfilled 



m. 



by a Langevin equation involving frictional and fiuctuation forces [2|, |3|. Another equivalent 
approach is the Fokker-Planck equation [4]. The latter constitutes a probabilistic description 
dealing, mainly, with evolved distribution functions. Here we should like to mention that 
the above approaehes have been suceessfuUy apphed to the fission process by eons.dering 
the deformation dynamics of atomic nuclei as a kind of Brownian motion |5i| . 

There exist, on the other hand, strong experimental evidences revealing that at low 
energies, the quanta, effeets are important at least .n the early stage of the reaofon H as 
is the case for the fast charge equilibration in heavy-ion collisions [71]. The description of 
dissipation phenomena at the quantum level is undoubtedly of fundamental significance, in 
particular at low temperatures where the quantum nature of the physical processes play a 
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crucial role |8| . It turned out that the direct quantization of the classical Langevin equation is 
quite problematic. Indeed, the usual quantization procedure can be carried out provided the 
equations of motion result from the classical Hamilton's principle. For this to be the case for 
the Langevin equation, the Lagrangian (Hamiltonian) should be explicitly time dependent. 
A first attempt in this direction consisted in the introduction of a time-dependent mass, 
which reproduces, in the classical limit, the frictional force |9[. It has been quickly recognized 
that this method exhibits fundamental difficulties. 

Later a nonlinear Shrodinger equation has been proposed by Kostin in order to describe 



dissipation in quantum systems. The above equation reads 10 1 
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where H is the Hamiltonian of the free system (i.e without dissipation), 7 is the friction co- 
efficient, and (.) denotes the expectation value. The origin of this equation can be elegantly 
explained using the Madelung hydrodynamical interpretation of Shrodinger equation: The 
gradient of the phase amplitude, given by the logarithm of the ratio of the wave function to 
its complex conjugate, is inversely proportional to the irrotational velocity field of the prob- 
ability fiuid. In addition to the nonlinearity of equation ([1]), which violates the superposition 
principle, its generalization to several degrees of freedom is not straightforward. 

It turns out that the natural way of describing dissipation phenomena in quantum me- 
chanics consists in dealing with the system as coupled to a large heat bath such that irre- 
versible energy fiows, from the former to the latter, take place |3|, |8|, [ll|, |l2| . Hence, the usual 
unitary group description of the evolution of quantum systems is not suitable for dealing 
with such processes since the irreversibility effectively introduces a preferable direction of the 
time. The latter fact can be accounted for using the notion of dynamical semi-groups, which 
are the generalization of unitary groups to non-Hamiltonian systems. This is, in particular, 
the case when the system of interest is opened, that is, when it is coupled to external systems 
or regarded as a part of a larger system. Recall that for a trace-preserving one-parameter 
semi-group, the group condition is merely replaced by: ^t+s 
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evolution equation for the reduced density matrix in the Schrodinger picture becomes 
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where H is the Hamiltonian of the system, and "D is a superoperator usually called the 
dissipator. Here we should like to mention that the coupling between the system and the heat 
reservoir is generally assumed weak so that the Markovian approximation is applicable 8|. 
This approximation is often used in quantum optics JjJ, and in studying nuclear fission 2|. 
The most general form of the dissipator I) giving rise to a quantum mechanical Markovian 
master equation was found by Lindblad [ij, [l5| ; it is given by 



nm) = ^ E([^^/^W' ^/] + i^^'PW^/])' 
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where the Ve are called Lindblad's operators; they depend on the system's variables and 
they model the effect of the environment. The dual evolution operator to ([2]) describes the 
time development of any Heisenberg operator A, namely, 

^ = l[H,Ait)] + ^5^([F/[i(t),f,] + [V,\Ait)m. (4) 

Depending on the system of interest, one can obtain axiomatically the master equation 
describing its damping by properly choosing the form of the Lindblad's operators. This 



axiomatic approach has 



3een extensively applied to the one- dimensional harmonic oscillator 



because of its simplicity 16|-|25| . The resulting master equation involves diffusion coefficients 
in coordinate and momentum which are of quantum nature |26| because of the fundamen- 
tal constraints they have to satisfy. The-two dimensional case has been considered in |27|, 
and [28|. It is, however, generally assumed that the off-diagonal elements of the quantum 
mechanical diffusion matrix are negligible. This, certainly, cannot be the case when the har- 
monic oscillators are strongly interacting with each other. In this paper we apply Lindblad's 
axiomatic approach to explicitly derive the elements of the diffusion matrix, for the gen- 
eral case of A^ harmonic oscillators, by imposing the condition that the asymptotic state is a 
Gibbs state. We should like to stress that the system we are considering here is of significant 
practical importance in nuclear physics, since the nuclear potential can be approximated, at 
least locally, by harmonic oscillators. In section |TT] we present a detailed derivation of the 
elements of the diffusion matrix and we investigate the constraints that should be satisfied 
by the transport coefficients. Section IIIII is devoted to the equations of motion for the ex- 
pectation values and variances of the coordinates an momenta operators. In Section IIVI we 
apply our results to the description of fast charge equilibration in deep-inelastic collisions, 
and we investigate the effect of dissipation on tunneling in sub-barrier processes. We end 
the paper with a summary. 

II. QUANTUM MECHANICAL TRANSPORT COEFFICIENTS FOR N COU- 
PLED HARMONIC OSCILLATORS 

The aim of this section is the derivation of the quantum mechanical diffusion coefficients in 
(mixed) coordinates and momenta for A^ coupled quantum harmonic oscillators by imposing 
a Gibbs asymptotic state for the evolved reduced density matrix. 



A. Explicit derivation 

The most general form of a A^-dimensional quadratic Hamiltonian in coordinates and 
momenta is given by 

N 
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+ ^ ^{i^kjqkQj + f^kjPkPj) + ^ fj-kjPkqj, (5) 



fe=i '' 



k^j k^j 

where the operators qt and pj satisfy the usual canonical commutation relations: 

[%,Pj\ = ihSkj, [qk,qj] = [Pk,Pj] = 0. (6) 

To ensure that the above Hamiltonian is physical the coupling strengths fikj, I'kj and k^j 
should satisfy certain conditions, some of which will be discussed later in particular cases. 
Here we just mention that, according to Onsager principle |29i], the last two coefficients 
should be symmetrical, that is: 

^kj = ^jk} f^kj = f^jk- \i) 

We require that the system relaxes to a steady state corresponding to A^ independent 
quantum harmonic oscillators in thermal equilibrium at temperature ksT = 1//3. This 
corresponds to complete thermalization of the system. Explicitly we have 

Peq = exp(-/3i^eq)/^, (8) 

where 

*^-E(^4"-2^2). (9) 

For the sake of generality, we assume that m^ and Uk are different from, respectively, M^ 
and Qk- The partition function, Z, can be easily calculated by rewriting i/eq in terms of the 
creation and annihilation operators, and then by taking the trace in the occupation number 

space; one finds that 

^ 1 

fc=l 



For convenience, we further introduce the following operators: 

^k = -j{QkPk+Pkqk), (12) 

Using the canonical commutation relations ([6]), we can show that the above operators are 
the generators of an 5*0(2, 1) Lie group, namely, 

[fl ff] = -2thwkfiSkj, [fl ff] = 2zhukf^Sk„ 

[flfj]=2ihoj,fl5,,. (14) 

Now, the Hamiltionian H and the density matrix peq may be expressed in terms of the 
operators T| as 

N , N N 



^ = Yl \^^^k + ^kfl + fkfl) + - ^{vkjQkQj + i^kjVkVj) + J] ^^kjVkqj (15) 

k=l ki^i k^j 

N 

K = Y.^3, (16) 
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Since peq = exp(— /3 ^^ Tl)/Z is a solution of the master equation ([2]) we should observe 
the following equality: 

C[Hs\ -Hs + ^ 5^(£[V^.]V/ - \C[V}V,] - ^V.^Ve) = 0, (18) 

where we have introduced the superoperator C whose action is defined by 



C[A] = exp[l3j2Ti)Aexp(^-^J2^i 

k k 



Taking into account the commutation relations (fT^ . a straightforward apphcation of the 
above equation yields: 

Clff] = cosh{2h(3ujj)fl + i smh(2h(3ujj)ff, (20) 

C[ff] = cosh{2H/3ujj)ff - i smh{2h/3ujj)f^. (21) 

The Lindblads's operators Ve may be expressed, in our case, as linear combinations of 
the coordinates and momenta operators. This is, in some sense, the quantum analogue of 
Hooke's law in classical mechanics. Consequently, we can write: 

TV N 

V, = Y.{a% + h%), V^^ = Y^{afqj + bfpj), £ = i;2iV, (22) 

j j 

where a^j and b^j are complex numbers. 

By making use of the canonical commutation relations ([6]) , together with equations ( !T9|) 
and f l20|) -f l2T|) . it can be shown that 



C[Hs]-Hs = 5^{[(54cosh(2n/3u;,)-l]-«/fcsinh(2n/3a;fc]f,i 
k 

+ [fk[cosh{2h/3ujk) - 1] + i5k smh{2h(3uJk)]f^ 
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QkQj 



+ - y^ \i^kjicosh.{h/3uk) cosh.(hf]ujj) — 1) — KkjmkmjUkOJj 

X sinh(^/3a;fc) sinh(^/3ci;j) + ifikjfnk<^k smh{h(3uk) cosh.{hf3u 

+ -/ \i^kjicosh.{hl3uk) cosh{h(3uj) — 1) 

2 ^^ L nikmjUkijJj 

X sinh(/i/3co'fc) sinh(/i/3cdj) + z ^ cosh(^/3a;fc) sinh(fi/3c(;j) p^pj 

rrikUk -I 

+ /^ \ fJ'kiicosh.(h(3uk) cosh.( hPuj) — 1) — ^-^ 

k^j ^''^'^ 

X smh.{h/3uJk) sinh(/i/3co'j) + 2i ( — smh.{hl3uk) cosh.{h/3ujj) 

\mkUJk 

+ Kkj-mjUj cosh.{hl3uk) smh.{hl3uj) j pkQj (23) 

The other terms of equation (TT8|) are easily calculated (see the appendix). In fact, the latter 
equation implies that the coefficients of all the involved operators should be equal to zero. 



As a result, by introducing the following notations: 



^.... = ^ReJ^^S'' D,,,^ = ^^^J2bi% (24) 
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D,,,^ = -^ReJ24*b'„ Xkj = -ImJ24*b'^, (25) 

e e 

akj = -ImJ^akO'p Vkj = -ImJ^bk^, (26) 

e e 

we obtain three independent sets of linear algebraic equations. The first one is as follows: 



[cosh(2/i/3a;fc) - 1] = -[Dq^q^mkUk ^^^^ ) sinh(;i/3wfc) cosh(/i/3wfc) 

sinh(/i/3a;fc) + h\kk\co'&\\{hfiuJk) + l); (27) 






smh{2hl3ujk) = --Dg^gJ cosh(/l/3a;fc) - 1) —smh{hl3uk) 



2mkUJk ^ ' rrikUJk 

+ ^^ sinh^ {h(3uk), (28) 

2hiikkmkUkSmh{2h/3uk) = Dj,^pJcosh{h/3uJk)-l) + hXkkmkUJkSmh{h/3uk) 

- Dq^qk^lulsinh^ihPuk), (29) 

4(cosh(2n/3a;fc)-l) = ^^^(l - cosh(n/3a;fc)') sinh(n/3a;fc), (30) 

The second one reads: 

nijUj 2mkDq^q.ujk ( ~1 + cosh.{h(5uJk) ) ( — 1 + cosh.{h(5ujj] 

+ ( hXjk + h^jk cosh.{h(5ujj) 1 smh.{h(5uJk) + ( hXkjmkCUk + h^kj'mkUJk cosh.{h/3ujk) 



-2Dp^pj smh.{hl3uk)) smh.{hl3ujj) = 0, (31) 

SDp^p. sinh^f — -^j sinh^f — -^j + rrijhujAXjk - iJ,jkCosh.{hl3uJk)] smh.{hl3uj) 
+mkhoJk smh{hl3uk) ( Xkj — fJ-kj cosh{h/3uj) — 2mjDg^q^ujj smh{hl3ujj) ) = 0, (32) 
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'' ^ ' cosh.{hj3LUj) — 1 ) smh.{h/3uk) — 2mjUjjDq^q^ ( cosh(^/3a;fe) — 1 ) sinh(fi/9ci;j) 



+fi< \ikj — Xkj cosh{h/3ujj) + cosh{hl3LUk)[Xkj — ^kj cosh(/i/3wj)] 
HjkrrijUj smh.{h/3uj) smh.{hf3uk 



. 0, (33) 

rUkUJk 



Finally, the equations of the third set can be written as: 

-2D,,,^ (cosh(n/3a;.) - l) (cosh(;^/3c.,) - l) + <^^^^^^^|^^^^^^ 
X smh.{h[5ujk) — rrijhuj ( akj + K^kj cosh.{hl3uk) ) smh.{hl3cuj) — 2 ( — ^-^ j 



\mkUJk' 
xDq.p^ smh{h(3uk) smh.{h(3uj) = 0, (34) 



-2Dq.p^ f cos\i{h(5ujk) ^ 1 ) ( cosh(/i/3a;j) — 1 j —nikhwk ( —auj + i^kj cos\i{h[5bjj) j 



X sinh(/l/3wfc) + ( j hrjkj + hvkj cosh.{hj3ujk) — '^f^kOOkD^^p. smh.{h(3ujk) 

X smh{hPuj) = 0, (35) 



hukj — hrjkj cosh.{hl3uj) + rrijUj ( —2Dg^p^ + hKkjmkUk smh.{hl3uk) ) smh{hl3uj) 
— cosh.{h(3uk) —hrjkj + hvkj cosh.{h[3uj) — 2Dq-p^mjUj smh.{h(3uj 
+2Dq^^p.mkUJk{cosh{hl3uj) - l) smh{hl3uk) = (36) 

From equations (1211) and ([2SD, it is clear that the quantities Dq^q^, ^PkPjy ^qwj ^^'^ ^kj are 
the multidimensional extension of the quantum mechanical diffusion and friction coefficients 
corresponding to the one-dimensional harmonic oscillator 
algebraic equations for the diffusion coefficients yields: 



Solving the above sets of 



r^ hi Xkk — l^kk \ , , hfiuk , „. 

Dq^Qk = o coth^-— , (37) 

2 y rUkUk J 2 

DpkPk = -rrikUJkihk + l^kk) coth , (38) 

Dp.,. = Dq,p, = --UkSk coth ^- = - (— - - 3^ jcoth ^-, (39) 



h ( Xjk — i^jk , , hf3cjk Xkj — fj'kj ,1 hf3(jjj 

yq^q - l^qq^ - - -^ '— COth \ '- ^^^ 

' ^'^ 4 \ nikUk 2 rrijUj 



Dq.q, = Dq^q^ = - ' ^/ COth ^— + ' ^^ COth -—^ , (40) 



T^ ^ h I .^ , _ h(3uk / , \ , h(3uj \ , ^ . , 

Dp^p^ = Dp^p^ = - I [Xjk + jijk)mkUk coth — h [Xkj + jXkjjrrijUj coth — --^ I , (41) 

^ ^ h ( rjkj + Vkj h(3uk h[5ujj\ , , 

D<i.p, = Dp^q^ = - coth — — + [akj - KkjjrrijUJj coth — — , (42) 

4 \ mkUJk / / 

which constitute the main result of this work. We can clearly see that the diffusion coeffi- 
cients (H0l) -( l42|) can by no means be neglected when the coupling constants Ukj and Kkj are 



comparable with, respectively, mkul and l/m^ (i.e. strong coupling). We shall discuss later 
the influence of these constants on the evolution in time of the relevant physical quantities. 
The diffusion coefficients form a N x N symmetrical matrix which we call the quantum 
mechanical diffusion matrix; from here on we shall denote it by -D. The diagonal elements 
of the latter along with the elements Dg^pi^ , corresponding to each degree of freedom, have 
the same form as those associated with the one-dimensional harmonic oscillator. This is 
a direct result of the quadratic form of the Hamiltonian H. Hence, in our model, the 
mutual interactions between the degrees of freedom do not affect the diagonal elements of 
the diffusion matrix. In particular, the fluctuation- dissipation theorem holds for every degree 
of freedom. Indeed, in the high-temperature limit, ksT ^ hujk, after a Taylor expansion of 
the trigonometric function, we obtain Einstein's relation: 

DpkPk = hhrnkksT, (43) 

where we have introduced the renormalized friction coefficient A^a; = ^kk + f'kk- 

Furthermore it is quite interesting to notice that all the other diffusion coefficients in 
mixed coordinates and momenta (140|) - (H2|) . are the arithmetic mean of two terms having 
the form of one of the coefficients fl37j) - fl39l) with, obviously, appropriate choice of both the 
phenomenological constants and the coupling strengths. Conversely, the latter coefficients 
may be obtained from the most general ones (H0 l) -( H2l) . by setting Uk = ojj, fijk = fikj = fJ'kk, 
Vkj = Mkfil, Kkj = 1/Mfc, and by observing that akk = Vkk = 0. Notice, also, that when 
Mk = nik, and Qk = ^k, then Dq^^p^, = 0, a value which is usually used in nuclear physics. 
By analogy to equation fH3|) . we may write, in the high temperature limit, 

Dp,p^=Kk,[^^^^)kBT, (44) 

where the friction coefficient in this case is given by 

1 



Kkj = {Xjk + fijk)mk + (Xkj + t^kj)mj . (45) 

Equation (l44l) is nothing but Einstein's relation for a fictitious one-dimensional harmonic 
oscillator whose mass parameter is equal to (m^ + mj)/2. When m^ = mj, jijk = fikj = 0, 
then Akj = {Xjk + Xkj)/2, i.e, the arithmetic mean of the friction coefficients Xkj and Xjk- 

It is worth mentioning that in the linear response theory, the one-dimensional diffusion 
coefficient Dpp is given in terms of the response function x {t) oi the operator that ensures 

10 



the coupling of the degree of freedom to the heat bath by pj 

/h6uj\ f°° 
Dpp{u) = cothf j / dtix {t)sm{ujt). (46) 

Relation ( 14T|) suggests that the multidimensional version of the above equation would be of 
the form 



^PkPj 2 



coth(^^-^j / dtix'jkit) sm{ukt) + coth{^-^j dtixlj{t) sm{u jt) , (47) 



where Xjk i^ ^^^ response function corresponding to the one- dimensional harmonic oscillator 
with frequency Uk and mass rrik which describes its coupling to both the heat reservoir and 
the oscillator with frequency Uj and mass rrij. 

B. Constraints on the values of the transport coefficients 

Let us first begin with briefly analyzing the conditions that should be satisfied by the 
coupling constants Ukj and Hkj appearing in the expression of the Hamiltonian operator 
H. In the special case where N = 2, and fikj = fJ'kk = 0, one can deduce, from simple 
mathematical considerations, that the kinetic coupling strength satisfies the inequality 



|ki2| < W . (48) 

V mi?7i2 

The coupling constant u^, on the other hand, is such that 

1 1/12 1 < ^/mrrrl^UlU2. (49) 

Consider now the case A^ = 3 with rris = m2 7^ mi, k,i2 = ^13 7^ ^23, ^12 = ^13 7^ ^23 and 
002 = 1^3 ¥" ^1- Then we should observe the following conditions: 



I I / ^ I I / n- + m2K23 / 1 , . 

«;23 < , k^l2 <a/^^ <\ , (50) 

m2 V 2mim2 V niim2 



\v2z\ < m2Ujj = mscol, |z/i2| < \ -miuf{u23 + m2(^l) < y/mrrn^uJiU2. (51) 

The quantum character of the diffusion coefficients we have derived above may be per- 
ceived from the fundamental constraints they have to satisfy. Indeed, taking into account 
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Cauchy-Schwartz inequality, we can infer from formulas (12^ -( 125|) that 

DqkqkDp.p, - Dg^p^ > —Xj,j, (52) 

(53) 

(54) 
(55) 

These conditions ensure the non-negativity of the density matrix at any moment of time. 
There exist in the literature, however, other sets of diffusion coefficients which violate these 



^qkqk^qjqj 


2 ^' 2 


^PkPk^PjPj 


Q lb C\ 

" ^PkP, - ~^^kj^ 


Oikk = Vkk = 


= 0. 



constraints. This is the reason for which they are usually called classical coefficients [26 1 
since a violation of the uncertainty relation may be observed at least at short times of the 
dynamics. 

In our multidimensional model, the friction coefficients cannot be freely chosen, in con- 
trast to the one-dimensional case where the friction coefficient is dealt with as a free param- 
eter which can be varied to reproduce the experimental data. As an illustration, consider 
the low temperature limit with yU^fc = /x^j = 0; then it is a matter of algebra to show that 



where 



\ I / L ^ UJ ^' / 



Ui - 2 



Tlkj + ^kj 



^mjmkOJkOJj 
Furthermore, in case where \kj = Xjk = 0, then 



+ y/mkm~uJkuJk{akj - Hkj) 



(57) 



\akj\ < \ '^^^ — , \Vkj\ < \/\kk>^jjmkmjUJkOJj, k y^ j. (58) 

y rrikmjUJkUJj 

The above conditions will be taken into account later in the numerical calculations. 



III. EQUATIONS OF MOTION 

In what follows we shall be interested in the evolution in time of the mean values and 
variances of the coordinates and momenta operators. The latter may be calculated once 
the density matrix p{t) is known. However, it is more convenient to work in the Heisenberg 
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picture. Using formula (jl]), it can be shown that the evolution in time of any Heisenberg 
operator F is given by 

kj I 



(59) 



X {F, qkQj} - ^kFqj ) + {2Dp^q^ + ihXjk) {F, {pj, Qk}} - 2{pjFqk + QkFpj 



- 2hXjk [ {F, pjqk} - 2qkFpj 



J / f' 



where {F, G} denotes the anticommuatator of the operators F and G. We recall here that 
the expectation values and variances are explicitly defined as 



or^(t) = tr(pF(t)), 
apcit) = ^trfp{F(t),G(t)})-ai.(t)aG(t) 



where tr(F) denotes the trace of the operator F. 
Let 



(60) 
(61) 



V(t) = Wq, (t) , (Tp, (t) , (Tg2 (t) , ap2 (t) , , ag^_^ (t) , ap^_^ (t) , a^^ (t) , a^^ (t) }' 



(62) 



and 



(^qiqiif) CTgiPlW (^qiq'zit) ^^qiPzi^) (^qiqwi^) ^giPiV W \ 



ait) 



'-^piqiV') '-^pipiV') '-^piqzV') '-^pipzV^) 



^qNqiV') ^qNPiv^) ^qNqzV^) '^qNPzV') 
y^PNqiV') ^pnpiV') ^PNq'iV') ^pnP2\'^) 
Then, by virtue of equation fl59|) . one can show that 

dV{t) 



dt 

da it) 

dt 



a, 



piqN 



(t) O-piPivW 



^qNqN V') ^qNPN V-) 
PNqNV^j ^pnpnv')/ 



a. 



MV{t), 

Mait) + a{t)M^ + 2D, 



(63) 



(64) 
(65) 
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where 



/ 



M 



-All +M11 


1 

ml 


— A12 + ^12 


-ai2 + K12 


— mi^i 


-All -/^ii 


'7l2 - 1^12 


— A21 — ^21 


-A21 +M21 


ai2 + Ki2 


— A22 + ^^22 


1 

m2 


-1112 - 1^12 


— A12 — /112 


-7712^2 


— A22 — /i22 



-AiTV + f^lN 
VlN - VlN 

~'^2N + fJ-2N 
V2N — V2N 



-aiN + Kin 
-^Ni — ^J■Nl 

-Oi2N + K2N 
-^N2 — fJ'N2 



-^Nl + ^J'Nl OiiN + Kin — A7V2 + fJ'N2 ■ ■ ■ 

\ —VlN — I^IN —^IN — I^-IN —V2N — V2N ' ' ' 

and D is the diffusion matrix: 

' ^qigi ^qipi ^qiq2 ^qip2 ^q\qN ^qiPN ' 

-L'piqi -^pipi ^piq-i ^PlP2 ^PiqN ^PlPN 

^qNqi ^qNPi ^qNqi ^qNP2 ^qnqN ^qnPN 

\^PNqi ^pnpi ^PNq2 ^pnP2 ^PNqN ^pnpn / 

The solution of equation flMl) is simply 

V(t) = exp(Mt)V(0), 



Aattv + ^J.NN :;k^ 

— niNiON ^^NN — ^-NN I 



m 



D 



(67) 



(68) 



whereas that corresponding to equation ( 165!) may be obtained by noting the following prop 
erty: 



dti 



*5e^*}=Ae^*Se^* + e^*5e^*C, 



(69) 



where the matrices A, B and C do not depend on time. Consequently, the time development 
of the matrix a is given by 



a{t) = exp(Mt)(a(0) - a) exp{Mtf + a, 



where the matrix a satisfies 



Ma + dM' +2D = 



(70) 



(71) 



In general the latter equation yields a set of —{2N + 1) = N{2N + 1) algebraic linear 
equations, the unknowns of which are the elements of the matrix a. 
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In our investigation we are assuming that the asymptotic sate of the system is the Gibbs 
state ([8]). Under this condition it can be verified that the expectation values and variances 
tend to 

hm ag^{t) = hm (Tp,(t) = hm ffp,g,(t) = 0, (72) 

lim (Tg^gd^) = (Tg^g.ioo) = - coth^-^, (73) 

t->-oo ZmkUJk ^ 

hm ap,p,(t) = (7p,p,(oo) = -mfcWfcCoth— -^, (74) 

I— >oo / / 

hm aq^q^it) = hm ap,p.(t) = hm (Tp^g^{t) = 0, k j^ j. (75) 

Hence it is possible to link the diffusion coefficients to the asymptotic variances by simple 
expressions. We have, for instance, 

1 



^IkPj ~ ^PjQk ~ 2 



(%j + z^fci)ffgfe<7fe(oo) + (afcj -fi:fej-)^P,P,(oo) • (76) 



Some remarks are in order here. First of all, the fact that the asymptotic expectation 
values (Tgj.(oo) and ap^{oo) are zero implies that the real part of all the eigenvalues of the 
matrix M should be negative, a fact that is equivalent to the condition exp(Mt) — )• as 
t — )■ oo. This, actually, imposes further conditions on the relevant parameters of the model. 
In particular we find that the matrix a is nothing but the asymptotic variance matrix, that 
is, a = cr(oo). 

Notice also that the generalized Heisenberg uncertainty relation 

itK,P,it)-^g,P,{t?>^ (77) 



a. 



qkQk 



4 



should be observed, since the operators p{t) and q(t) satisfy the usual canonical commutation 
relation at any moment of the time. It has been shown, however, that when the fundamental 
constraints imposed on the diffusion coefficients are not satisfied, then it may happen that 
the inequality (!77|) is violated at certain interval of the time. 

The expectation value of the Hamiltonian H can be calculated using equation f l59|) . The 
resulting formula is quite cumbersome, and we shall not display it here. Nevertheless, by 
direct calculation one can verify that 

E = lim ti{pH{t)) = -yuk coth ^^. (78) 

k 
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IV. APPLICATION TO HEAVY-ION COLLISIONS 

In what follows, we shall apply the results obtained above to the description of the motion 
of a dinuclear system (DNS) in the charge and mass asymmetry coordinates 

Z1-Z2 iVi-iV2 ,_., 

Here Zi, A^^i and Z2, N2 are, respectively, the charge number and neutron number of the 
nuclei. The advantages of the DNS come into play in the description of the various reaction 
channels in heavy ions collisions, such as fission and fusion of atomic nuclei. 

Based on the work of Hahn et a/ (7|, Sandulescu et al 30J, the authors of 3l| proposed 



an analytically solvable quantum-mechanical model describing the charge and mass distri- 
bution in heavy-ion collisions. There, the investigation consists in solving the following 
time-depending Schrodinger equation: 



h^ d h^ d 1, 2 1 2 



2Mzzdri'z 2MNNdrij, 2 " '^ 2 



ip{Vz,VN,t) 



d 
ih—ip{riz,VN,t), (80) 



where Mzz and Mj^jy are mass parameters, kz and k^- are stiffness parameters, and kzN 
is the coupling constant. They have, however, neglected dissipation by assuming that the 
mass and charge asymmetry degrees of freedom are isolated from the other collective and 
intrinsic degrees of freedom. They also assumed that the neutron and proton mobilities are 
uncorrelated; this is the reason why there is no momentum-momentum coupling in the above 
equation. Before we proceed further, note that the quadratic form of the Hamiltonian in (I8OI) . 
valid only for nearly grazing collisions, was obtained from an expansion around the energy 
surface minimum 1]^ = rjz = oi the potential energy between two colliding nuclei, which is 
defined as the sum of the usual liquid-drop energy, the Coulomb contribution due to charged 
protons, the rotational energy (proportional to the square of the total angular momentum), 
and the proximity nuclear potential. The last two contributions depend strongly on the 
relative distance between the two interaction partners. 

Our task here is to investigate the effect of the transport coefficients on the dynamics of 
the compound nuclear system by introducing through Lindblad's formalism damping effects 
. Attention will be given to the influence of the coupling between the collective degrees of 
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freedom on the evolution in time of the expectation values and variances (see reference 27 1 
for further discussion). 



(a) 




-2.0 - \ 



Q.O 






(b) 



t(10 --S) 



t(lO-^-s) 



FIG. 1: (Color online) Time evolution of: (a) azzit), and (b) azpj^{t) for zero off-diagonal elements 
of D (dashed lines) and nonzero off-diagonal elements of D (solid lines). The dot-dashed line 
represents the asymptotic value corresponding to the Gibbs state. The parameters are Mzz = 
Mnn = 461. 6344/1 VMeV, tiuoz = 2.9468MeV, hwN = 2.9288 MeV, vzN = -1869 MeV, hXzz = 
HXnn = 2 MeV, T = 5 MeV, azz{0) = W-\ a^^pM = ^V(4fTzz(0)), a^N = IQ-^, ap^p^(O) = 
ft^/(4(TArAf(0)); all other parameters are set to zero. 



A. Illustrative calculations 



Unless otherwise stated, the model parameters we shall use in the sequel are those 
corresponding to the system ^^^Xe + ^^^Sn. The stiffness parameters kz and kj^ are 
found to be equal to 4009 MeV and 3960 MeV, respectively 3l|. The average mass pa- 
rameters, calculated within the framework of the hydrodynamical theories, are given by 
Mzz = M]siN ~ 461.6344/i^/MeV. This corresponds to a value of the angular frequencies 
of hujz ~ 2.9468 MeV, hujN ~ 2.9288 MeV. The calculation gives a value of 3739 MeV for 
the coupling constant kzN, which implies that vzn = ^nz = —1869 MeV. It can easily 
be checked that these values satisfy the condition (149 p . The friction coefficient A^^ (with 
k = N, Z) has the dimension of the angular frequency Uk, they are, in general, of the same 
order. More precisely, due to fast charge equilibration, we should have 2AArAr > un and 
2A^^ > Uz- The initial value of the mass and charge asymmetries can easily be calculated; 
one can find that riz{0) = 0.0385, and //at = 0.0067. 
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FIG. 2: (Color online) Time evolution of the expectation values aj\f{t), (Jz{t), (7p^{t) and crp^{t) 
for different values of the coupling constant unz] ^nz = 3000 MeV (solid line), vnz = —1968 MeV 
(dashed line) and vnz = -3000 MeV (dot-dashed line). Here T = 2 MeV, ap^(O) = ap^{^) = 0; 
the other parameters are the same as figured) 

In order to illustrate the importance of the off-diagonal elements of the diffusion matrix 
we display in figure [1] the evolution in time of the variances (Tzzif) and o-zpf^if) for both 
zero and nonzero off-diagonal elements. We see that the asymptotic values of the variances 
when the off-diagonal elements are set to zero do not correspond to the Gibbs state ([H]) 
as indicated by the dot-dashed line in the above figure. Therefore, we conclude that the 
behaviour of the dynamics of the compound nuclear system depend strongly on the values 
of the diffusion coefficients. 

An example of the development in time of the expectation values and variances of the 
charge and neutron asymmetry coordinates for different values of the coupling strength v^z 
is displayed in figures |2] and El It can easily be seen that the expectation values decay faster 
for large positive values of the latter parameter. The decay gets slower as we decrease vnz to 
negative values, which implies that the evolution of the mean values depends on the sign of 
the coupling constant. This fact manifests itself even though the system is nearly symmetric. 
The differences in the evolution in time of the proton and neutron asymmetry coordinates 
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FIG. 3: (Color online) Evolution in time of the variance cr^zit) for different values of the coupling 
constant vnz', ^nz = 3000 MeV (solid line), vnz = —1968 MeV (dashed line) and vnz = —3000 
MeV (dot-dashed line). Here T = 2 MeV; the other parameters are the same as figure [TJ 

may be explains by the fact that the dynamics is sensitive to the initial expectation values. 
Notice that the situation is, however, slightly different regarding the evolution in time of 
the variances cnz and (Jpj^p^r (their values quantify the correlation between the two degrees 
of freedom). In this case the curves corresponding to couphng constants having the same 
magnitude but with different signs are symmetrical, see figure [3l The investigation shows 
that the above result does not hold when the initial values of a^z and cTpj^p^ are different 
from zero. The other variances do not change much when varying unz- 

Let us now investigate the effect of the friction coefficients and the kinetic coupling 
constant knz on the behaviour of the DNS. Figure H] displays the time dependence of the 
expectation values for different values of Xnn and Xzz- As expected we see that the decay 
of the above quantities is less appreciable for small values of the friction coefficients. This 
result does not qualitatively change with nonzero values of the remaining model parameters. 
It is also found that except ajyz (see figure [5]), and CTpr^pz, the other variances are not much 
affected by changing the values of the friction coefficients. The time development of the 
centroids for different values of knz is illustrated in figure El Once again we find that the 
decay is faster for large values of the coupling constant, whereas the correlation between 
the proton and neutron asymmetry coordinates becomes larger, as shown in figure [71 When 
the initial values of a^z and Opj^p^ are zero then the above quantities are symmetrical with 
respect to the change of the sign of ktvz (see figure [3] for a similar situation). The other 
variances are robust with regard to the variation of the latter constant. 
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FIG. 4: (Color online) Time evolution of the expectation values afq{i\ o-z{t), (Tp^{t) and crp^{t) 
for different values of the friction coefficients Xnn and Xzz] fiXjyN = hXzz = 3 MeV (solid line), 
HXnn = hXzz = 2 MeV (dot-dashed line) and hX^N = hXzz = 1-6 MeV (dashed line). Here 
T = 2 MeV, h'zN = —1869 MeV, crp^{0) = crpj^{0) = 0; the other parameters are the same as 
figure [Ij 

B. Comparison with experimental data 

Now we are going to assess the results of our model by comparing them with the exper- 
imental data obtained by Schiill et al [32|. For this reason we have calculated the ratio of 
neutron to proton variances and the correlation coefficient Xnzi defined by 

(yNz{t) 



XNz{t) 



(81) 



\/crNN{t)(Jzz{t) 

for the reaction ^^^Xe + ^^^Sn. The results are displayed in figure [HI The theoretical curves 
(solid lines) were obtained for a value of the friction coefficients Xnn = ^zz = 2 MeV/h, 
with azN = —otNZ = 33 x 10^*^ MeV~^ s~^ which gives the best fit to the experimental data 
without violating the fundamental constraints on the transport coefficients. We can see that 
with the diffusion coefficients ( P7|) - (H2]) . the theoretical values of the ratio crj^]\r(t)/(Tzz(t) are 
in quite good agreement with the experimental ones in the interval < t < 10^^^ s. On the 
other hand, though not in perfect agreement with the experimental outcomes, the model 
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FIG. 5: (Color online) Evolution in time of the variance cr]\fz{t) for different values of the friction 
coefficients Xnn and Xzz', f^^NN = fiXzz = 3 MeV (solid line), fiX]\f]y = HXzz = 2 MeV (dot- 
dashed line) and hXjyN = hXzz = 1-6 MeV (dashed line). Here T = 2 MeV; the other parameters 
are the same as figure [H 



gives improved results for the correlation coefficient as compared with those of 3l|] and 33 |. 
This difference may be explained by the nonzero value of the parameter azN, which is 
responsible for the creation of momentum-momentum correlations between the charge and 
mass asymmetry coordinates, even though the motion of the neutrons and protons was 
assumed to be initially uncorrected {unz = 0). We have further checked the validity of 
the above results by coupling the neutron and proton asymmetry coordinates to the relative 
motion of the nuclei. It turns out that the only difference between the two cases resides 
in a small diminution of the value of the friction coefficients. It is worth mentioning that 
the ratio of neutron to proton variances does not depend much on the value of a^z- Also, 
nonzero initial values of aj^z do not significantly improve XNz{t)- 



C. Penetration enhancement due to dissipation in sub-barrier processes 

As a second application, in connection with the description of the fusion process, we 
now investigate the penetration of Gaussian wave packets through a potential barrier ap- 
proximated by a two-dimensional inverse harmonic oscillator. It should be stressed that 
there exist no metastable states for this kind of potentials; a more appropriate one would 
be composed of a potential well smoothly linked to a parabolic barrier. 

The solutions of the equations of motion for this case may be obtained by simply making 
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FIG. 6: (Color online) Time evolution of the expectation values iTAr(i), crz{t), (Tp^{t) and crp^{t) 
for different values of the coupling constant knz', knz = 33 x 10^^ MeV~^ s^^ (solid line), Kjyz = 
20 X lO^s MeV"^ s-2 (dot-dashed line) and knz = (dashed line). Here T = 2 MeV, CTp^(O) = 
(Tpj^(O) = 0; the other parameters are the same as figured) 

the replacement 

Uk -^ ioJk (82) 

in the expression of the matrix M [see equation ( 1^ ]. In the subsequent calculations we 
use the diffusion coefficients ( l37|) - fH2|) . In general, the density matrix in coordinates space 
is given by the integral of the Wigner distribution function of the system with respect to 
momentum variables 3J|, namely. 



p{<ii,<i2,---<iN;t) 



OO /"OO 



OO J — OO 



N 



Yl dpkW{qi,pi, ga, P2 ■ ■ ■ Qn, Pn] t). 



(83) 



fc=i 



It has been shown using path integral techniques that for the harmonic oscillator, if the 
initial density matrix is Gaussian, then it remains Gaussian at any moment of the time 23] . 
This means that the Wigner function can be determined by a simple substitution of the 
time-dependent variances and expectation values, namely. 



N 



W{qi,pi,q2,P2---qN,PN;t) = (27r) 2 A/det(cr ^{t))exp 



-iZ 
2^ 



V{t)fa-\t){Z-V{t)) 

(84) 
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FIG. 7: (Color online) Evolution in time of the variance crNzit) for different values of the coupling 
constant knz; k-nz = 33 x 10^^ MeV~i s"^ (^^^y^,^ Yme), knz = 20 x 10^® MeV"^ s~2 (dot-dashed 
line) and knz = (dashed line). Here T = 2 MeV, (7p^{0) = crpjy(O) = 0; the other parameters are 
the same as figured! 

with Z = {qi,Pi,q2,P2- ■ ■ (1n,Pn}^, and c^^it) denotes the inverse of crit). In figure [9] 
we display the time development of the probability density function fl83|) for the case of a 
two-dimensional parabolic barrier. One can see that at short times p{qi,q2,t) gets sharper 
as compared with its initial shape; obviously, this is accompanied by an increase of its 
height since its area should be constant, (more precisely, the integral of the above quantity 
over the whole space should be equal to unity.) The height of the distribution function 
decreases with time whereas p{qi,q2,t) spreads in coordinates space, to become centered 
around (crg^(oo), (7^2(00)) at sufficiently long times. 

The probability of finding the packet to the right of the barrier in the gi direction is given 
by (note that the top of the barrier is located at the origin) 



Pit) 



dq2 / dqip{qi,q2]t) 




00 00 00 

dP2 / dpi / dq2 / dqiW{qi,pi,q2,P2;t). (85) 



P{t) is used here to quantify the penetrability through the parabolic barrier. Notice that 
in one dimension, enhancement of the tunneling was found for large values of the friction 
coefficient 23|]. Here we shall investigate the effect of dissipation when other degrees of 
freedom are considered. 

In figure [TOl the penetration probability is shown as a function of time for different values 
of the friction coefficient A22- It can bee seen that in the interval < t < 20 x 10^^^ s, P{t) is 



23 



N 






-i 




0.6 








- 


0.4 


: ] 










^______- 


N 

^ 0.2 










L-^ 


: 


0.0 


--^ 


'- 


— -- 




-0.2 


"1 


- 



0.0 



0.5 



1.0 

t (10"^^s) 



1.5 



0.0 



0.5 



1.0 
t (10"^^s) 



1.5 



FIG. 8: (Color online) Experimental data [32] (dots) along with the theoretical curves (solid 
lines) corresponding to the proton to neutron variances ratio (left) and the correlation coeffi- 
cient XAfz(i) (right) as functions of time for the reaction ""^^^Xe + ""^^"^Sn. Here the parameters are: 
Mzz = Mnn = 461.6344nVMeV, fvjJz = 2.9468 MeV, hwN = 2.9288 MeV, vzn = -1869 MeV, 
hXzz = nXNN = 2 MeV, T = 0.02 MeV, azz{0) = lO'^, ap^p^(O) = nV(4azz(0)), aNN = lO'^, 
'^Pivpjv(O) — ^^/(4o"ArAf(0)), (7j\[z{0) = 0, azN = 33 X 10^^ MeV~^ s~^; all other parameters are set 
to zero. 

almost the same for all values of A22- Then the curves spread apart from each other, to tend 
to certain asymptotic values which depend, in turn, on the dissipation rate. Indeed, we see 
that the greater the value of the friction coefficient, the larger the asymptotic penetrability, 
as clearly indicated in figure [TOl Thus the tunneling through the barrier in the gi direction 
is enhanced by the dissipation in the other degree of freedom. Quite surprisingly, we find 
that for sufficiently large values of A22, the wave packet is trapped near the top of the barrier 
(Pit) ^ 0.5). 

We have also studied the dependence of the tunneling on the coupling constant uu- It 
turns out that P{t) is inversely proportional to the latter parameter. We also found that 
the penetrability increases with the temperature which can be explained by the increase of 
the values of the diffusion coefficients. 

V. SUMMARY 

In this paper we have axiomatically derived the multidimensional diffusion coefficients 
for a set of N coupled harmonic oscillators using Lindblad's approach. The only assumption 
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FIG. 9: (Color online) The probability density function p{qi,q2',t) at different values of time. The 
parameters are mi = 2.bh^/yieV, m2 = GO^i^/MeV, hwi = 1.7 MeV, nw2 = 0.6 MeV, 1/12 = 7 
MeV, ^Aii = 2.5 MeV, h\22 = 0.6 MeV, T = 0.1 MeV, ^^^^^(O) = 0.4, ap^p^{<d) = ^V(4cTgigi(0)), 
(0) = 7 X 10-2, ap^pM = ^V(4^g2<?2(0)), ^91 (0) = -6, (jg,(0) = 0, a^M = ^^, ^^^2(0) = 0; 



'^92 '72 

all other parameters are set to zero. These values correspond to an initial total energy Eq = 
-115.9 MeV. 

we have made is the existence of a Gibbs steady state for the system under consideration. 
It turns out that the general form of the off-diagonal elements of the diffusion matrix is 
given by the arithmetic mean of two terms, each having the form of one of the diffusion 
coefficients corresponding to the one- dimensional harmonic oscillator. Furthermore, we have 
shown that the fluctuation-dissipation theorem holds for both diagonal and off-diagonal 
coefficients in momentum coordinates. The elements of the friction tensor are found to be 
not independent. This is due to the fundamental constraints on the diffusion coefficients. 
We have derived the equations of motion for the expectation values and variances, and 
solved them for arbitrary values of the coupling strengths, without having recourse to any 
perturbative treatment. This is, indeed, one of the advantages of the investigated model. 
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FIG. 10: (Color online) Probability of penetration through the barrier P{t) as a function of time 
for different values of the friction coefficient A22- The parameters are mi = 2.5/i'^/MeV, m2 = 
GOfi^/MeV, fvjji = 1.7 MeV, ^2 = 0.6 MeV, 1/12 = 7 MeV, ^An = 2.5 MeV, T = 0.1 MeV, 

o'gigil'JJ = 0-4, o-p 

—6, O"g2(0) = 0, o"pj(0) = 9/i, crp2(0) = 0; all other parameters are set to zero. These values 

correspond to an initial total energy Eq = —115.9 MeV. 



.,,,,(0) = 0.4, cTp,p,(0) = h'/{Aa,,,M), ^<?2.2(0) = 7 x 10"^ a^^^^M = h'/{Aa,,,M). ^<?i(0) 



We have applied our results to the description of mass and charge asymmetry coordinates in 
deep-inelastic collisions. We find that the expectation values of the coordinates and momenta 
do not depend on the temperature, in contrast to the variances which are temperature 
dependent. The decay of these quantities is faster for large values of both the friction 
coefficients, and the coupling constants. The correlation between the degrees of freedom 
is more appreciable when the coupling is strong. This was confirmed by comparing the 
theoretical results with the experimental data. It is also shown that dissipation in one degree 
of freedom enhance the tunneling in sub-barrier processes. In conclusion, the model is quite 
interesting, in the sense that it is exactly solvable; further extensions and investigations may 
be carried out. 



26 



Appendix A 



The result of applying the superoperator £ [see equation (fT9|) ] to the operator V^Vi is 
riven by 

£[V^^Ve] = V J [|4|2 cosh\hf3uk) —{4% + aibi*) cosh{hf3uk) smh{hPuk) 

^—^ 1 L TTLhUJh 



,2, ,2 



1 r / W'? \ 

sinh^f^/^cjfc) y\, + af 6^ cosh^(^/3a;fc) + i ( la^pmfcCUfc ^-^^ ) 

J L V muuju) 



X cosh ( h(3ujk ) sinh {h(3uk) + a^ 6^* sinh^ ( /i/^cufc 



Pkqk + 



IJL\2 

■rrikUJk' 

I.U* „„„T,2/ 



46f cosh''(;i/3c<;fc) 



+ i[\au\ rrikUk 



j cosh(/i/3a;fc) smh{hl3uk) + a-khk sinh {h/Suk 



+ 






PkQk 



|6^| cosh {h/3uJk) + inikOOkicik^k + '^fc*^fc) cosh(/i/3co'fc) sinh(/i/9a;fc) 



- |4pm^w^sinh2(n/3a;fc)]g2l+^J J4*a^^ cosh(;i/3 



Wfe) cosh(/i/3ci;j) 



— i cosh.{hl3uJk) smh.{hl3uj j] 



nikmjUkUJj 



bk< 

i — coshi hf3LL!j) smh.(h(3LL!k) 

rUkUk 



smh.{h/3uk) smh.{hf3uj 



PkPj + 



bl^b'j cosh.{hPuk) cosh.{hl3ujj) 



+ ibka^jTrijUj cosh.{hl3u)k) smh{h(3uj) + ia^kb''jmkUJk cosh{h(3uj) smh.{hl3uk) 



— a^a,mkmjijjkijjj sm\i{h(3ijJk) sm\i{fi(3ijj j) 



qkQj + 



(4*65 + aib'^ cosh(n/3a;,) 



X cos}i{hj3uJk) + i{a^k^] + al.af')7njUj cosh{hl3uk) smh.{hl3uj) — 



rrikUJk 



ifnj^j ,,i^ I ,ei*\ 



X {biW^ + bibf) smh{hPuk) cosh{h^Uj) + ^^^{bl*a] + blaf) sui\i{n(5ujk 



VflkUk 



X smh.{hl3ujj] 



PkQj 
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